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Oh| Abstract. It is considered that the effective interaction between any two quarks in 

^ . a baryon can be approximately described by a simple harmonic potential. Also, it 

' is made use of the nonrelativistic approximation since the effective (constituent) 

1/^ , masses of quarks are not small. The problem is firstly solved in Cartesian coor- 
dinates in order to find the energy levels irrespective of their angular momenta 
and then it is also solved in polar cilindrical coordinates for taking into account 

^ , the angular momenta of the levels. By making a comparison between the two 

' solutions the energies and the corresponding angular momenta (and parity) of al- 

I most all baryon levels are described. The agreement with the experimental data 

. is quite impressive. The solution in Cartesian coordinates also produces some 

' very important figures for the sizes of baryons and for the harmonic oscillator 

, constant which is clearly related to confinement. 



Introduction 



As is well known there are several important works that deal with the calcu- 
lation of the energy levels of baryons. One of the most important ones is the 
(~| , pioneering work of Gasiorowicz and Rosner [1] which has calculation of baryon 

levels and magnetic moments of baryons using approximate wavefuncions. An- 
other important work is that of Isgur and Karl [2] which strongly suggests that 
' non-relativistic quantum mechanics can be used in the calculation of baryon 

^ , spectra. Other very important attempts towards the understanding of baryon 

spectra are the works of Capstick and Isgur [3], Bhaduri et al.[4], Murthy et 
al.[5], Murthy et al.[6], and Stassat et al.[7]. An important work attempting to 
describe baryon spectra is the recent work of Hosaka, Toki and Takayama [8] 
published in 1998. This last work arrives at an important equation which had 
already been deduced by De Souza a long time ago, in 1992 [9]. Other works 
by De Souza published before 1998 include it [10,11]. 

The effective potential between any two quarks is not known and because 
of this several different potentials are found in the literature. In particular the 
harmonic approximation using a harmonic central potential has been widely 
used. Since the three quarks of a baryon are always in a plane it is assumed 
that the effective potential between any two quarks of a baryon can be given by a 
linear harmonic oscillator. The motion of the plane of quarks is not considered 
in this paper. This is a calculation quite different from those found in the 
literature. 
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2 Calculation in Cartesian coordinates 



In the initial calculation we use normal cartesian coordinates which, of course, 
does not consider the angular momentum of the system, that is, it does not take 
into account the symmetries of the system. But this section is very important 
because it calculates the energy levels. In the next section we will link each level 
to its angular momentum. Considering the work of Isgur and Karl [2] as to the 
use of non-relativistic quantum mechanics and using a linear harmonic oscillator 
potential [8,9] we can write the Hamiltonian in normal cartesian coordinates as 

where we have used the fact that the three quarks are always in a plane. The 
above equation may be resolved into a sum of 6 equations 



-^(Ei- -c.,Ci V' = 0, (2) 



which is the equation of a single harmonic oscillator of potential energy 50;,^^ 
and unitary mass with E = ^i- 

The general solution is a superposition of 6 harmonic motions in the 6 normal 
coordinates. The eigenfunctions V'i(Ci) are the ordinary harmonic oscillator 
eigenfuntions 

mi) = A^..e-(«*/2)«?//,,(^e,), (3) 

where Ny. is a normalization constant, at = vi/Ti and Hy.{^/ai^ij ^ Hermite 
polynomial of the u^th degree. For large ^^ the eigenfunctions are governed by 
the exponential functions which make the eigenfunctions go to zero very fast. 
The energy of each harmonic oscillator is 

E^ = hl^i{v^ + i), (4) 

where Vi = 0,1,2,3, ... and i/j is the classical oscillation frequency of the normal 
"vibration" i, and Vi is the "vibrational" quantum number. The total energy of 
the system can assume only the values 

E{vi,V2,V3, ...ve) = hi'i{v-i + ^) + hi'2{v2 + ^) + ■■■hveive + ^). (5) 

As was said above the three quarks in a baryon must always be in a plane. 
Therefore, each quark is composed of two oscillators and so we may rearrange 
the energy expression as 

E{n,m,k) = hviin -\- 1) + hv2{m + I) + hv^ik -i- 1), (6) 

where n = ui + ?;2, m = 173 + 774, = W5 + ug. Of course, n, m, k can assume the 
values, 0,1,2,3,... We may find the constants from the ground states of some 
baryons. They are the known quark constituent masses taken as to„ = ma = 
0.31Gev, rus = O.SGev, nic = 1.7Gev,m6 = 5Gev and mt = 174GeV. 
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The states obtained with the above Hamihonian arc degenerate with respect 
to isospin so that our calculation does not distinguish between nucleonic and A 
states, or between E and A states. In the tables below the experimental values 
of baryon masses were taken from reference 12. 

Let us start the calculation with the states ddu(neutron), uud(proton) and 
ddd(A~), uuu (A++) and their resonances. All the energies below are given in 
Gev. Because m„ = m^, we have that the energies calculated by the formula 

Er,,m,k = 0-Sl{n + m + k + 3) (7) 

correspond to many energy states. The calculated values are displayed in Table 
1. The last column on the right is a rough classification which will be cleared 
up in the next section. 

The energies of the particles A and S, which are composed of uus and uds 
are given by 

i^n,m,fc = 0.31(n + m + 2) + 0.5(fc + l). (8) 

The results are displayed in Table 2. The agreement with the experimental val- 
ues is excellent. For the E°{uss) and E~{dss) baryons the energies are expressed 

by 

i^n.m.fe = 0.31(n+l)+0.5(m + fc + 2). (9) 

See Table 3 to check the agreement with the experimental data. In this case 
the last column is almost empty due to a lack of experimental data. 
In the same way the energies of Cl{sss) are obtained by 

-£^n,m,fe = 0.5(n + m + A; + 3). (10) 

The energies are displayed in Table 4. The discrepancies are higher, of the 
order of 10% and decreases as the energy increases. This is a tendency which is 
also observed for the other particles. This may mean that, at the bottom, the 
potential is less flat than the potential of a harmonic oscillator. 

The energies of the charmed baryons(C = +1) A+, E+ and E° are 

given by 

En,m,k = 0.31(n + m + 2) + 1.7{k + 1). (11) 

The levels arc shown in Table 5. 

For the charmed baryons (C = +1) 5+ and Sj! we have 

En^m,k = 0.31(n + 1) + 0.5(m + 1) + 1.7{k + 1). (12) 

The results are displayed in Table 6. 
As for the fl^, its energies are 

i^n,m,fe = 0.5(n + m + 2) + 1.7(fc + l). (13) 

Table 7 shows the results of the energy levels. 

In all tables below Ec is the calculated value, Em is the measured or exper- 
imental value and the Error is 

Em — Ec 



Error 



X 100%. 



Ec 

As is easily seen in each table many levels may be predicted. 
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n, m, k 


Ec{Gev) 


-EM(Gev) 


Error(%) 


L2I.2J 


Parity 


0,0,0 


0.93 


0.938(7V) 


0.9 


Pll 


+ 


n -\- TTi k = 1 


1.24 


1.232(A) 


0.6 




+ 


n -\- m -\- k = 2 


1.55 


1.52(A^) 


1.9 


Di3 




n -\- m -\- k = 2 


1.55 


1.535(iV) 


1.0 


i ± 




n -\- m -\- k = 2 


1.55 


1 6fA) 


3.1 


-^33 




n -\- m -\- k = 2 


1.55 


1.62(A) 


4.5 






77 -1- 777 -4- h — ^ 

11/ 1 III/ 1 11/ KJ 


1.86 


1.90(A^) 


2.2 


^ io 




f) _|_ irn -1-/^—3 

9 V \ 9 9 V \ 9\j t_J 


1.86 


1.90(A) 


2.2 


*-^3i 




Tl -\- Tfl -\- k = ^ 


1.86 


1.905(A) 


2.4 


-^35 




Tl -\- TTI -\- k — 2f 


1.86 


1.91(A) 


2.7 






71 + 772 + A) = 3 


1.86 


1.92(A) 


3.2 


P33 




Ti -\- vn -\- k = A. 


2.17 


2.{]8(N) 


4.1 


^ io 




fl -\- TTl -\- k = 4: 


2.17 


2M{N) 


3.7 


5*11 




Tl -\- TTI -{- k = 4: 


2.17 


2.10{N) 


3.2 


Pn 




n -\- m -\- k = A 


2.17 


2.15(A) 


0.9 


S'si 




77 + 777 + /c = 4 


2.17 


2.19(A^) 


0.9 






77 + 777 + A; = 4 


2.17 


2.20(A^) 


1.4 


Di5 




77 + 777 + fc = 4 


2.17 


2.20(A) 


1.4 


G37 


_ 


n + m + k — 4 


2.17 


2.22(7V) 


2.3 


Hig 


+ 


n + m + k = 4 


2.17 


2.225(iV) 


5.5 


G19 


- 


n + m + k = 5 


2.48 


2.39(A) 


3.6 


F37 


+ 


n + m + k = 5 


2.48 


2.40(A) 


3.2 


G39 




n + m + k = 5 


2.48 


2.42(A) 


2.4 


Hs.n 


+ 


n + m + k = 6 


2.79 


2.7(iV) 


3.2 


^1,13 


+ 


n + m + k = 6 


2.79 


2.75(A) 


1.4 


-^3,13 




n + m + k = 1 


3.10 


3.100(Af) 





-^1,15 


7 


n + m + k = 8 


3.21 


? 


? 


? 


7 


Tl + m + k = 9 


3.72 


? 


? 


? 


7 


n + m + k = 9 


4.03 


? 


7 


7 


7 















Table 1. Baryon states N and A. The energies Ec wc;rc! (;aleulated ac- 
cording to the formula En,m,k = 0.31(77 + 777 + fc + 3) in which n, m, k are 
integers. Em is the measured energy. The error means the absolute value 
of {Ec — Em)/Ec- We are able, of course, to predict the energies of many 
other resonances. 
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State(n, m, k) 


EciGev) 


Em {Gey) 


Error(%) 


L2I,2J 


Parity 


0,0,0 


1.12 


1.116(A) 


0.4 


Ui 


+ 


0,0,0 


1.12 


1.193(1;) 


6.5 


Pii 


+ 


71 + 771 = 1 , k=0 


1.43 


1.385(1!) 


3.2 


lo 


+ 


n -\- m = 1, k=0 


1.43 


1.405(A) 


1.7 






71 + 771 = 1, k=0 


1.43 


1.48(i;) 


3.5 


? 


? 


0,0,1 


1.62 


1.52(A) 


6.2 






0,0,1 


1.62 


1.56(1!) 


3.7 


? 


? 


1 


1.62 


1.58(i;) 


2.5 






1 


1.62 


1.60(A) 


1.2 


Pm 




1 


1.62 


i.62(i;) 









0,0,1 


1.62 


i.66(i;) 


2.5 




+ 


oAi 


1.62 


1.67(S) 


3.1 






0,0,1 


1.62 


1.67(A) 


3.1 






77 + 777 = 2, k=0 


1.74 


1.69(A) 


2.9 






77 + 777 = 2, k=0 


1.74 


1.69(S) 


2.9 


7 


? 


n + m = 2, k=0 


1.74 


1.75(E) 


0.6 


Sii 




77 + 777 = 2, k=0 


1.74 


1.77(S) 


1.7 


Pll 


+ 


77 + 777 = 2, k=0 


1.74 


1.775(S) 


2.0 


Di5 




n + m = 2, k=0 


1.74 


1.80(A) 


3.4 


Soi 




77 + 777 = 2, k=0 


1.74 


1.81(A) 


4.0 


Poi 


+ 


77 + 777, = 2, k=0 


1.74 


1.82(A) 


4.6 


-Fo5 


+ 


71 + II) = 2. k=() 


1.74 


1.83(A) 


5.2 







Continues on next page 



5 



State(n, m, k) 


Ec{Gev) 


-EM(Gev) 


Error(%) 


L2I,2J 


Parity 


n -\- m = 1 , k=l 


1.93 


1.84(S) 


4.7 


io 


+ 


n -\- m = 1, k=l 


1.93 


1.88(11) 


2.6 


Pn 




n -\- m = 1 , k=l 


1.93 


1.89(A) 


2.1 




+ 


n + m = 1, k=l 


1.93 


1.915(S) 


0.8 




+ 


n + m = 1, k=l 


1.93 


1.94(S) 


0.5 


Di3 




n + TO = 3, k=0 


2.05 


2.00(A) 


2.5 


? 




n + TO = 3; k=0 


2.05 


2.00(1)) 


2.4 


Si 1 




n + TO = 3, k=0 


2.05 


2.02(A) 


1.5 




+ 


n + TO = 3, k=0 


2.05 


2.03(S) 


1.0 


Fn 


+ 


n + TO = 3, k=0 


2.05 


2.07(S) 


1.0 


Fl5 


+ 


n + TO = 3, k=0 


2.05 


2.08(S) 


1.5 


Pl3 


+ 


0,0,2 


2.12 


2.10(S) 


0.9 


Gn 


- 


0,0,2 


2.12 


2.10(A) 


0.9 


Go7 


- 


0,0,2 


2.12 


2.11(A) 


0.5 




+ 


n + TO = 2, k=l 


2.24 


2.25(S) 


0.5 


? 


? 


n + TO = 4, k=0 


2.36 


2.325(A) 


1.5 


^03 


- 


n + TO = 4, k=0 


2.36 


2.35(A) 


0.4 


Ho9 


+ 


n + TO = 1, k=2 


2.43 


2.455 


2.5 


? 




n + TO = 3, k=l 


2.55 


2.585(A) 


1.4 


? 


? 


0,0,3 


2.62 


2.62(E) 





? 


? 


n + TO = 5, k=0 


2.67 


to be found 


? 


? 




n + TO = 2, k=2 


2.74 


to be found 


? 


? 




n + TO = 4, k=l 


2.86 


to be found 


? 


? 




n + TO = 1, k=3 


2.93 


to be found 


? 


? 




n + TO = 6, k=0 


2.98 


3.00(E) 


0.7 


? 


? 


n + TO = 3, k=:2 


3.05 


to be found 


? 


? 




n = m = 0, k=4 


3.12 


to be found 


? 


? 




n + TO = 5, k=l 


3.17 


3.17(E) 





? 


? 


n + TO = 2, k=3 


3.24 


to be found 


? 


? 

















Table 2. Baryon states E and A. The energies Ec were calculated according 
to the formula m.k = 

0.31(n + TO + 2) + 0.5(A: + 1). Em is the measured 
energy. The error means the absolute value of {Ec — Em)/Ec- 
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State(n, m, k) 


-Ec(Gev) 




Error(%) 


L2I,2J 


Parity 


0,0,0 


1.31 


1.315 


0.5 


Pll 


+ 


1,0,0 


1.62 


1.53 


5.6 


Pl3 


+ 


1,0,0 


1.62 


1.62 





? 


? 


1,0,0 


1.62 


1.69 


4.3 


? 


? 


n=0, m + k = 1 


1.81 


1.82 


0.6 






2,0,0 


1.93 


1.95 


1.0 


? 


? 


n=l, m + k = 1 


2.12 


2.03 


4.2 


? 


? 


n=l, m + k = 1 


2.12 


2.12 





? 


? 


n=3, m = k = 


2.24 


2.25 


0.5 


? 


? 


n=0, m + k = 2 


2.31 


2.37 


2.6 


? 


? 


n=2, m + k = 1 


2.43 


to be found 


? 


? 


? 


n=4, m = k = 


2.55 


2.5 


2.0 


? 


? 


n=l, m + k = 2 


2.62 


to be found 


? 


? 


? 



Table 3. Baryon states S. The energies i?c were calculated according to 
the formula En^m,k = 0.31(n + 1) + 0.5(m + A; + 2). Em is the measured 
energy. The error means the absolute value of {Ec — Em)/Ec- The state 
!B(1530)Pi3 appears to be the lowest state of the composite S |+) tt. Its decay 
is in fact Htt. 



State(n, m, k) 


^c(Gcv) 


£'m(Gcv) 


Error(%) 


0,0,0 


1.5 


1.672 


11.7 


n + m + fc = 1 


2.0 


2.25 


12.5 


n + m + k = 2 


2.5 


2.47 


1.2 


n + m + A; = 3 


3.0 


to be found 


? 











Table 4. Baryon states il. The energies Ec were calculated according to 
the formula En,m,k = 0.5(n + m + k + 3), and Em is the measured energy. 
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State(n, m, k) 


-Ec(Gev) 


Em (Gey) 


Error(%) 


0.0,0 


2.32 


2.285(Ac) 


1.5 


n + m = 1, k=0 


2.63 


2.594(Ac) 


0.1 


n + m= 1, k=0 


2.63 


2.627(Ac) 


0.01 


n + m = 2, k=0 


2.94 


to be found 


? 











Table 5. Baryon states Ac and Ec. The energies Ec were calculated ac- 
cording to the formula En,m,k = 0.31(n + to + 2) + l.l{k + 1). The state 
with energy 2.63 MeV had already been predicted in another version of this 
work. The experimental levels 2.594 MeV and 2.627 MeV have confirmed 
the theoretical values. It appears that the level Sc(2.455) is a composition 
of the level (0, 0, 0)(that is the 2.285 Ac) with a pion as is also inferred from 
its decay. 



State(n, to, A;) 


Ec{GeY) 


i^M(Gev) 


Error(%) 


0,0,0 


2.51 


2.47(S+ 


1.6 




2.51 


2.47(s;; 


1.6 


1,0,0 


2.82 


to be found 


? 


0,1,0 


3.01 


to be found 


? 











Table 6. Baryon states Sc. The energies Ec were calculated according to 
the formula En,m,k = 0.31(n+ 1) +0.5(m-|- 1) + 1.7(A;-|- 1). Em is the mea- 
sured energy. The recently found level Sc(2645) is probably a composition 
of the regular level S+ with a pion as its decay confirms. 
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State(n, m, k) 


-Ec(Gev) 


Em (Gey) 


Error(%) 


0.0,0 


2.7 


2.704(r2'j 





n + m = 1, k=0 


3.2 


to be found 


? 


n + m = 2, k—0 


3.7 


to be found 


? 











Table 7. Baryon states Oc- The energies Ec were calculated according to 

the formula En,m.k = 0.5(n + m + 2) + 1.7(/c + 1). The energy of the level 
(0, 0, 0) above shown had been predicted in other versions of this work. 



We may predict the energy levels of many other baryons given by the for- 
mulas: 

• ucc and dec, £^„,„j,fe = 0.31(n + 1) + 1.7(m + A: + 2); 

• sec, En,m,k = 0.5(n + 1) + 1.7(m + k + 2)- 

• ccc, En,Tn,k = l-7(n + m + fc + 3); 

• ccb, En,m,k = l-7(n + m + 2) + 5(fc + 1); 

• ebb, En,m,k = l-7(n + 1) + 5(m + k + 2); 

• ubb and dbb, En^rn,k = 0.31(n + 1) + 5(m + k + 2) ; 

• uub, udb and ddb, En,m,k = 0.31(n + m + 2) + b{k + 1); 

• bbb, En^m.k ~ 5(n + m + fc + 3); 

• usb and dsb, En^m,k = 0.31(n + 1) + 0.5(m + 1) + 5(A; + 1); 

• sbb, En,m,k = 0.5(n + 1) + 5(m + fc + 2); 

• scb, K.m.fc = 0.5(n + 1) + 1.7(m + 1) + 5(fc + 1); 

• ucb, En,m,k = 0.31(n + 1) + 1.7(m + 1) + 5(fc + 1); 

• ttt, En,m,k = (174 ± 17)(n + m + fc + 3); 

• and all combinations of t with u, d, c, s and b. 
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3 Calculation in polar cilindrical coordinates 

In order to address the angular momentum and parity we have to use spherical 
or polar coordinates. Since the three quarks of a baryon arc always in a plane 
we can use polar coordinates. Wc choose the Z axis perpendicular to this plane. 
Now the eigenfunctions are angular momentum eigenfunctions (of the orbital 
angular momentum). Thus, wc have three oscillators in a plane. Considering 
that they are independent the radial Schrodinger equation for the stationary 
states of each oscillator is given by [13] 



where is the quantum number associated to L^- Therefore, what we have 
is the following: three independent oscillators with orbital angular momenta 
Li, L2 and which have the Z components Lzi,Lz2 and L^s, in the plane 
containing the quarks. Of course, the system has a total orbital angular mo- 
mentum L = Lx + L2 + L3 and there is a quantum number li associated to 
each Li. The eigenvalues of the energy are given by-'^^ 

E = (2ri + |mi| + l)hPi + (2r2 + |m2| + l)hu2 + (2r3 + jmsl + l)hu3 (15) 

in which ri, r2, rs = 0, 1, 2, 3, ... and \mi\ = 0, 1, 2, 3 , li. Comparing the above 

equation with the equation 



we see that n = 2ri + m — 2r2 + \m2\, k = 2rs + Im^]. 

Let us recall that if we have three angular momenta Li, L2 and described 
by the quantum numbers li,l2, h the total orbital angular momentum L will be 
described by the quantum number I given by 



where h > |mi|, I2 > \m2\, I3 > |m3|. 

Taking into account spin we form the total angular momentum given by 
J = L + S and the quantum numbers of J arc j = I ± s where s is the spin 
quantum number. As we will see we will be able to describe almost all baryon 
levels. 

3.1 Baryons N and A 

Let us begin the calculation with the particles N and A. Wc will classify the 
levels by energy according to Tabic 1. The first state of N is the state (n = 
0, TO = 0, A; = 0) with energy 0.93 GeV. Therefore in this case h = I2 = I3 = 
and then / = 0. Hence this is the positive parity state Pn and we have 




(14) 



E{n, m, k) = hui{n + 1) + hu2{'m + 1) + hu3{k + 1) 



/l+/2 + /3>i>||«l-/2|-/3| 



(16) 
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I N A Parity 



0.938Pii ? + 



The second energy level (1.24 GeV) which is the first state of A has n + 
m + k = 1. This means that 2ri + \mi \ + 2ri + \mi\ + 2r3 + \m-i\ — 1. Thus, 
\fn\ \ + |to2| + Imal = 1 and + ^2 + ^3 > 1, and we can choose the sets \m\\ = 
1, \mi\ = Imal = 0; |toi| = Imal = 0, |m2| = 1; |mi| = 1, |m2| = Imal = 0, and 
= 2,^2 = ^3 = 0, or h = 2,/i = Za = 0, or still ^3 = 2,Zi = ^2 = which 
produce I = 2 and thus the Icvvl 



I N A Parity 



2 ? I.232P33 + 



In the third energy level (1.55 Gev) n + m + k = 2 = 2ri + |mi| + 2r2 + 
|m2| + 2r3 + |m3|. This moans that |mi| + |m2| + Imal =2,0 and we have the 
sets of possible values of Zi, ^2, ^3 





2,0,0 0,2,0 0,0,2 1,1,0 


1,0,1 


0,1,1 0,0,0 


I 


2 2 2 0,1,2 


0,1,2 


0,1,2 



in which the second column presents the values of I that satisfy the condition 
+ ^2 + ^3 > 2, 0. There are thus the following states 



I 


N 


A 


Parity 





1.44Pn, I.7IP11 


? 


+ 


1 


1.5355ii, 1.65511 
1.52£)i3, 1.70i?i3, 1.675£)i5 


1.62531 
I.7OD33 




2 


1.68Fi5,1.72Pi3 


I.6P33 


+ 



because we can have j = 1/2 = + 1/2 = 1 - 1/2; j = 3/2 = 1 + 1/2 = 
2 - 1/2; j = 5/2=1 + 3/2 = 2 + 1/2. 

The fourth energy level (1.86 Gev) has n + m + A; = 3 = 2ri + |mi| + 2r2 + 
|to2| + 2r3 + |m3| which makes |mi| + |m2| + |m3| = 3, 1 and h + 12 + h > 3, 1. 
We have therefore the possibilities 



h, h, h 


/ 


h, h, h 


I 


^1, h, h 


I 


3,0,0 


3 


2,0,1 


3,2,1 


0,2,1 


3,2,1 


0,3,0 


3 


1,0,2 


3,2,1 


1,0,0 


1 


0,0,3 


3 


1,2,0 


3,2,1 


0,1,0 


1 


2,1,0 


3,2,1 


0,1,2 


3,2,1 


0,0,1 


1 
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and the states 



I 


N 


A 


Parity 


1 


2.08L'i3 


1.90531 

I.94D33 




2 


1.90Pi3,2.00Fi5, 1.99i^i7 

2.00^35,1-95^37 


I.9IP31, I.92P33, 1.905^35, 
? 


+ 


3 


? 


1.93£>35 





In the fifth energy level (2.17 Gov) n + m + k = i = 2ri + |mi| +2r2 + |m2| + 
2r3 + |m3| which yields |mi| + |m2| + Ims] = 4, 2, and /i + /2 + ^3 > 4,2,0. We 
can then have h = I2 = I3 = {I = 0) and also 



^1, h, h 


I 


^1, I2, h 


I 


4,0,0 


4 


2,2,0 


4,3,2,1,0 


0,4,0 


4 


2,0,2 


4,3,2,1,0 


0,0,4 


4 


0,2,2 


4,3,2,1,0 


3,1,0 


4,3,2 


2,0,0 


2 


3,0,1 


4,3,2 


0,2,0 


2 


1,3,0 


4,3,2 


0,0,2 


2 


1,0,3 


4,3,2 


1,1,0 


2,1,0 


0,3,1 


4,3,2 


1,0,1 


2,1,0 


0,1,3 


4,3,2 


0,1,1 


2,1,0 



and hence the states 



I 


N 


A 


Parity 





2.IOP11 


? 


+ 


1 


2.085ii, 2.20i:>i5 


2.15531, 2.35D35 




2 


7 


2.39P37? 


+ 


3 


2.19Gi7, 2.25Gi9 


2.2OG37 




4 


2.227^19 


2.3if39 


+ 



In the sixth energy level (2.48 Gev) n + m + k = 5 = 2ri + |mi| + 2r2 + |m2| + 
2r3 + |m3| which produces \mi \ + \m2\ + Im^l = 5, 3, 1 and h + h + h > 5, 3, 1. 
We have then the possibilities 
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5,0,0 


5 


0,5,0 


5 


0,0,5 


5 


4,1,0 


5,4,3 


4,0,1 


5,4,3 



= h,l2,h 


/ 


- 1,4,0 


5,4,3 - 


- 1,0,4 


5,4,3 - 


- 0,4,1 


5,4,3 - 


- 0,1,4 


5,4,3 - 



^1, h, h 



3,0,0 


3 


0,3,0 


3 


0,0,3 


3 


3,1,0 


4,3,2 


3,0,1 


4,3,2 



- hjhth 


I 


- 1,3,0 


4,3,2 


- 1,0,3 


4,3,2 


- 0,3,1 


4,3,2 


- 0,1,3 


4,3,2 



Thus we identify the states 



I 


N 


A Parity 


2 


7 


2.39F37 + 


3 


? 


2.406-39 


4 


? 


2.42fl"34i + 


5 


2.60/1, 11 


? 



The seventh energy state (2.79 Gev) has n + m + k = 6 = 2r\ + \m\\ + 
2r2 + |m2| + 2r3 + |m3| which produces \m.i\ + \m.2\ + \m^\ = 6,4,2,0 and 
^i + ^2+^3>6,4,2,0. We have then the possibilities below 



h,h^ I?, 


0,0,0 0,G,0 


0,0,6 


5,1,0 


5,0,1 


1,5,0 


1,0,5 


I 


6 6 


6 


6,5,4 


6,5,4 


6,5,4 


6,5,4 



and the states 



I 


N 


A 


Parity 


4 


? 


? 


+ 


5 


? 


2.75/3,13 




6 


2.7/^1,13 


2.95/^:3,15 


+ 



3.2 Baryons S and A 

Now let us do the calculation for S and A. According to Table 2 the first 
energy state (1.12 Gev) is (n = 0, m = 0, fc = 0) and hence we can have 
h = 0, Z2 = 0, Z3 = which yields I = and the states 



I s 


A 


Parity 


1.193Pii 


I.II6P01 


+ 
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In the second energy level (1.43 Gev) n + m = l,k = which makes 2ri + 
|mi| + 2r2 + \m2\ = 1 and 2r3 + Imaj = 0. This actually makes \mi \ + \m2 \ = 1 
and Imsl = 0. That is, we have the condition h + h > ^, h > which allows us 
to choose the possibilities 



ll,l2,h 1,1,0 


1,0,1 


0,1,1 


I 2,1,0 


2,1,0 


2,1,0 



that produce the states 



I 


E 


A 


Parity 





1.385Pi3 


? 


+ 


1 


? 


1. 405501 




2 




7 


+ 



and the state 1.48S is either S13, Sii{l = 1) or Fi5{l = 2). 

In the third energy level (1.62 Gev) n = m = 0,A; = 1 and we have |mi| = 
0,1^2! = and Imal = 1. That is, we have the condition li > 0,^2 > 0, 
^3 > 1 which allows us to choose h = I2 = OJs = = I3 = IJs = 0; 
Zi = 0, Z2 = ^3 = 1, and the states 







A 


Parity 





I.66P11 


I.6OP01 


+ 


1 


1.625ii 
1.58£)i3 


1.675oi 
1.52£>o3 




2 


? 


? 


+ 



and then the state 1.56S is probably Fi5( I = 2). 

The fourth energy level (1.74 Gev) has n + m = 2 = 2ri + \mi \ + 2r2 + \m2\ 
and k = 2rs + jmsl = 0, and thus we obtain \mi \ + |to2| =2,0 and \m3\ = 0. 
Hence we have the condition h + I2 > 2,0 and I3 > 0. We can then choose 
li = 2,l2 = h = 0;li =13 =0,^2 = 2;li = ^2 = 1,^3 = and thus we can 
identify the states 





E 


A 


Parity 





1.77Pn 


I.8IP01 


+ 


1 


1.75511 
1.67£'i3, 1.775£)i5 


l.SOSoi 
1. 692^03 




2 


? 


1.82Po5 


+ 
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and then the level 1.69S is probably Fi5(Z = 2). 

In the fifth energy level (1.93 Gev) n + m = 1 = 2r\ + \m\\ + 2r2 + |m2| 
and k = \ = + Ims], and thus we obtain \m\\ + \mi\ = 1 and \m^\ = 1. 
Hence we have the condition /i + ^2 > 1 and > 1. We can then have the sets 
l\ = 1, Z2 = 0, ^3 = 1; /i = 0, ?2 = 1, ^3 = 1- Both yield I = 2, 1, and we can 
identify the states 



I 


E 


A 


Parity 





1.84Pn,1.84Pi3 


1.89Po3 


+ 


1 


1.94_Di3 


1.83_Do5 




2 






+ 



The sixth energy level (2.05 GeV) has n + to = 3 = 2ri + \m\ \ + 2r2 + |to2| 
and fc = = 2r3 + Imsj, and thus wc obtain \m\ \ + |to2| = 3, 1 and |to3| = 0. 
Hence we have the condition + ^2 > 3, 1 and Iz > 0. Wc can then have the sets 
h = 2,l2 = 1, ^3 = 0; /i = 1, /2 = 2, ^3 = which make / = 3, 2, 1, for h + 12 > S 
and the sets h = 1, ^2 = 1, ^3 = 0; = 1, ^2 = 1, ^3 = which make / = 2, 1, 0, 
for Zi + ^2 > 1- We can identify the states 



I 


E 


A 


Parity 





2.O8P13 


? 


+ 


1 


2.005ii 


? 




2 


2.07Fi5,2.03Fi7 


2.02^07 


+ 


3 


? 


? 





In the seventh energy level (2.12 GeV) n = = 2ri + |toi|,to = = 
2r2 + \m2\,k = 2 = 2r3 + |to3| and thus |toi| — 0, |to2| = 0, |to3| = 2,0. Hence 
we have the condition ?i > 0, ^2 > and I3 > 2, 0. We can then choose the sets 
h =0,l2 = 0, ^3 = 2; /i = 0, I2 = OJs = 3 which make Z = 3, 2, and the states 



/ 


E 


A 


Parity 


I = 2 


? 


2.IIF05 


+ 


1 = 3 


2.IOG17 


2.IOG07 





Unfortunately, the angular momenta of the other energy levels have not been 
found but they can surely be explained according to what was developed above. 
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3.3 Bciryons S 

For these baryons only some angular momenta are known. The first energy level 
(1.31 GeV) has n = 0,m = O.l = which make = ^2 = '3 = and / = and 
is thus a P state. Therefore we obtain 



/ S Parity 



I.3I8P11 + 



In the second energy level (1.62 GeV) n = 1 = 2ri + |toi|,to = = 2r2 + 
|to2|,A: = = 2r3 + |to3| and thus \mi\ = 1, |to2| = 0, |to3| = 0. Hence we 
have the condition ?i > 1,^2 > and I3 > 0. We can then have the sets 
^1 = 1, Z2 = 0, 13 = 0; h = 1, Z2 = 1, ^3 = which make / = 2, 1, 0, and the states 



I 




Parity 





1.53Pi3 


+ 


1 ■! 


2 


? 


+ 



and thus the two levels 1.62 and 1.69 are probably either S, D or F states. 

The third energy level (l.SlGeV) has n = = 2ri + |toi|,to + fc = 1 = 
2r2 + \m2\ + 2r3 + |m3| and thus \mi\ = 0, |m2|+ Imaj = 1. Hence we have the 
condition li > 0,l2+ h > 1- We can then have the sets h =0,^2 = 1,^3 = 0; 
li = 0, Z2 = 0, ^3 = 1 which make / = 1, and the state 



I S Parity 



I = 1 1.82D13 



In the fourth energy level (1.93GeV) n = 2 = 2ri + |mi|,m = = 2r2 + 
|m2|, A: = = 2r3 + Imsl and thus |mi| = 2,0, |m2| = 0, jmsl = 0. Hence we 
have the condition li > 0,2,^2 > and ^3 > 0. We can then choose the set 
h = 2,l2 = 0, Z3 = which produces I = 2, and the state l.QSGeV is probably 
an F state. 

3.4 Relation between energy and angulcir momentum 

Prom Eqs. 15 and 20 we have 

E = (2ri + |mi| + l)hui + (2r2 + \m2\ + l)hu2 + (2r3 + Imgl + l)/iz/ 
Zi + ^2 + ^3 > ^ > 11^1 - ^2! - ^3] with li > |mi|,Z2 > |'^^2M3 > l^sl- 
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in which li, I2 and I3 are the quantum numbers of the angular momenta Li, L2, 
and L3, and mi,rn2, ms are the quantum numbers of their projections on the Z 
axis, respectively. Therefore, we clearly see that levels with large energies have 
large angular momenta as is quite evident from the experimental data. 



4 The sizes of baryons 

The solution in Cartesian coordinates is also useful for calculating in a quite 
simple manner the average size of a baryon. As is known the average potential 
energy of each oscillator is half of the total energy, that is, 

<lk^.'>-^iv.+ '^) (17) 

but since there are two directions for each quark in the plane there actually are 
two oscillators per quark and thus we have the potential energy Eg associated 
to each quark 

Eg = hugini + 1) (18) 

where rii = 0, 1,2,3, ... and hvq is the constituent quark mass constant. Thus 
taking into account Eq. 06 and the above fact on the relation between the total 
energy and the potential energy for an oscillator it can be written that 

E{n,m,k) = hui{n + 1) + hi^2{m + 1) + hi/sik + 1) = 
= 2 X (< ifcir?!^ > + < ifc2r72^ > + < ^hm^ >) = (19) 

= < kiTJi"^ > + < fc2??2^ > + < fc3»73^ > 

where = Cij^ + ^ifc^ which j and k are the two orthogonal directions of 
the two oscillators. One can then make the association 

hiy,{n + 1) = < kirji^ > (20) 

and hence the average radius 7?. of a baryon can given by 

/ , \l/3 

n{n,m,k) = [^/<r]l'^ ><r]2'^ ><r]3'^ >) = 

- (21) 

fei fe2 fea 

It is quite obvious that the application of the above formula should be first 
done to the proton. In the fundamental level n = m = k = and hvi = hv2 = 
hi/ 3 = O.SlGeV, and making the reasonable supposition that ki = k2 = k^ = k, 
thus 

7eo=7e(0,0,0) = y^. (22) 

If one uses for the average size of a proton the figure of •\/0.72fm = 0.85fm [14] one 
has k « 0.5GeV/fm^ which is a very reasonable figure because if it is multiplied 
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by the characteristic distance of Ifm (of course) the constant k' w 0.5GeV/fm 
is obtained which is quite close to the value of the constant K used in the QCD 
motivated potential [15], [16] 

Vqcd = -CF-^ + Kr (23) 

which is assumed to be of the order of IGeV/fm. 

From Table 1 one has that for n = m = fc = 2 the energy of a proton is about 
2.80GeV which gives an average radius of about 1.39fm and hence one sees that 
the size of a baryon docs not change much with the its energy. Therefore it 
can be said that the smallest radius of a proton is about O.Sfm and the largest 
radius is or the order of 1.4fm. 

For the ground states of E~ and S~ reference [15] gives, respectively, the 
radii •\/0.54fm = 0.73fm and yOTiSfm = 0.66fm. In terms of quarks S~ is dus 
and therefore one should have ^2 = fcs and ki k 0.5GeV/fm^ and 

1/3 



TZin, m,k)=l | hMn + l)><hMn^^+^l)><hus{k+l) ^^4) 



Using the above value it is obtained that 

/ /n.3i X n.3i X n.5\ 
0.73 = , 



0.5 X {ksf J 

which yields ks = 0.80GeV/fm^. Prom Table 2 it is seen that for n + rn = 5, 

fc = 1, the energy is 3.17GeV which is the highest energy level up to now. If one 
takes, for example, n = 2, m = 3 one has an average radius of about 1.24fm. 
Also it is found that the average radii of E~ are much smaller than those of 
the proton for levels with the same quantum numbers n, m, k. 

Now one can turn to S~ which in terms of quarks is dss. Then it is expected 
to have the same fcs « 0.80GeV/fm^ (two of them) above and a new fc, which 
can be called kgs- Using the ground state radius of (0.66fm) one obtains 
kss ~ 1.47GeV/fm^. For the excited states the average radius (in fm) is thus 

1/3 



«(n,"a-) ^ (^^ °.31(~+1)X0.5(^J^1)X0.5(. + 1) ^ 

which for the highest known excited state 2.55GeV (n = 4, to = A; = 0) gives 
7^(4,0,0) « 1.48fm. Using the value kss ~ 1.47GeV/fm2 the radius of the 
ground state of is estimated to be about 0.58fm. 

Putting together the above values the very important table below (Table 8) 
is obtained for the constant k (which is a sort of constant of confinement) in 
terms of the pairs of interacting quarks. 

The table shows that k increases with the reduced mass of the pair of interact- 
ing quarks. When the data are fitted to a polynomial up to second order in 
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the reduced mass of the pair of interacting quarks the foUowing polynomial is 
obtained 

kin) = 0.1188 - 1.7561m + 28.6508^^ (26) 

It is interesting that the coeficient of the last term is quite large and thus the first 
derivative increases very rapidly with /j,. As more massive quarks are considered 
the degree of the polynomial may increase but just to have a lower bound one can 
calculate the value of k for the interaction between two top quarks. The above 
formula gives k = 7416GeV/fm^. If the above data are fitted to a polynomial 
with a higher degree, for example, A:(/i) = A + Bji^ + C/x"^, the following values 
are obtained: A = -0.0268, B = 23.3794, and C = 0.2278. Since the value of C 
is small and B is of the same order of 28.6508, the first polynomial (Eq. 26) is 
a good approximation. If it is used for obtaining the k between quarks u and c 
one has k{uc) w 1.53GeV/fm2, k{sc) w 3.7GeV/fm2, k{cc) w 19GeV/fm2. And 
then one has that the radii of the ground states of the charmed baryons A+, 
S++, E+ and E° are about 




0.31 X 0.31 X 1.7 
0.5(1.53)2 




which is not so small due to the influence of the interaction between the two u 
quarks. As to Oc its ground state has a radius 

1/3 




, /0.5 X 0.5 X 1.7 , 
^-"'-^ 1.53(3.7)2 ) -'-''^ 

and the ground state of the ccc baryon has the quite small radius of just 




Rccc \/ TTTT ~ 0.3fm. 



Since the value of k{cc) « 19GeV/fm2 was obtained by means of an extrapola- 
tion the above figure of Rccc should be taken as a crude approximation. 

In the case of the ttt baryon an even cruder number is gotten for its ra- 
dius because its value for k is expected to be larger than the above figure of 
7416GeV/fm2, but it is instructive anyway to calculate its order of magnitude 
which in this case produces an upper bound for its radius. Therefore one can 
say that the radius of the ground state of the ttt system 
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Rttt < \l = 0.15fm. 

*" V 7416 

which is a very important number just because the top quark is the most massive 
quark. 

Since in this work the motion of the plane where quarks are sitting was not 
taken into account conclusions can not be drawn on the shape of baryons using 
the above figures. 
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5 On the spin-orbit interaction 



We clearly notice that the splittings of some levels are caused by the spin-orbit 
interaction. For example, consider the states 1. 905*31 and 1.91i-3i of A which 
differ by the values of / = 1 and I = 0, respectively. Since we are assuming a 
harmonic potential and as the spin-orbit term is proportional to we can 
approximately write 

AEsL ~ C < 5 . L >= C + 1) - 1(1 + 1) - sis + 1)] (27) 

for N and A baryons which have quarks with equal masses. Using for the above 
case j = 1/2, s = 1/2 we find C « 5MeV which shows that the influence of the 
spin-orbit interaction is small. Considering the levels I.9IP31 and I.92P33 we 
find C « 3.3MeV which is of the same order of the above C. The same holds 
in the case of the other baryons: for example, consider the states 1.755*11 and 
I.77P11 of S or the states 1.805oi and I.8IP01 of A. We see that there is a small 
energy difference between these states. 

6 Conclusion 

The simple model presented above which considers that a baryon is composed 
of three nonrelativistic quarks produces very importante results. First of all it 
describes quite well almost all energy levels of baryons with the proper assign- 
ment of the states as to parity and angular momentum. And the calculation 
also yields some reasonable figures for the sizes of baryons and for the important 
harmonic oscillator constant k which is directly related to confinement. 
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